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NOMENCLATURE 


Aj = coefficients of PadI approximant 

fR = wing aspect ratio 

Bi = coefficients of partial fraction of Pad6 approximant 

bp = local span at the section containing the transducer as 

shown in Fig, 17(b) 

= reference chord 
c* = semi chord length of the airfoil 

C(k) = Theodorsen function 

Cc(k) = generalized Theodorsen function 

C] = two-dimensional oscillatory lift coefficient 

= two-dimensional steady state lift coefficient 
Cm = oscillatory pitching moment due to sinusoidal gusts 

F = Fourier transform 

F(k) = C(k)[Jo(k)-iJ,(k)] + iJi(k) 

h(s) = impulse lift function due to nuclear blast 

Jo(k), Ji(k) = Bessel functions of the first kind 

k = w c / U , reduced frequency 

k' = k/0.61, effective reduced frequency 

L = Laplace transform 

L = oscillatory lift distribution 

L'(k) = in-phase component of oscillatory lift coefficient 

L"(k) = out-of-phase component of oscillatory lift coefficient 

M = freestream Mach number 

r = ik, Laplace transform variable 

s = Ut/c* or Ut/c, non-dimensional distance parameter 

S = reference wing area 


V 



t 

T(k) 

U 

V(k) + iW(lc) 
Wa 


wa* 

Wa 

Wg 

Wg* 

Wg 

Wo 

X 

>«le 


x,z 


a 

6 

Bt 

6{t) 

0 

P 

<P 

Hs) 

( 1 ) 


= time coordinate 

= L* (k) + iL"(k) , oscillatory lift coefficient due to 
sinusoidal gusts 
= frees tream velocity 
= l-T(k), [see Eq. (30)] 

= vertical velocity component of the vortex sheet simulat- 
ing the airfoil in gust field 
= amplitude of Wg 

= Wa*/U 

= vertical velocity component of gust 
= amplitude of vertical velocity component of gust 
= w^*/U 

= vertical velocity component of a sharp-edged gust 
= x/c* or x/c, non-dimensional x-coordinate 
= x-coordinate of wing leading edge 
= reference point for gust phase 

= distance penetrated into a step gust, [see Fig. 17(b)] 
= Cartesian coordinate system [see Fig. 1] 

= angle of attack 
= 

= root of the polynomial in the denomenator of Pad6 

approximant, which is always real 
= delta function 

= blast intercept angle (see Fig. 20) 

= frees tream density 

= phase difference (in degrees) of oscillatory lift 
= indicial lift due to step gusts 
= radian frequency associated with the gust wavelength 
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1. INTRODUCTION 


Estimation of the response of an aircraft due to atmospheric gusts 
has been the subject of numerous investigations from the viewpoint of 
producing useful data on the induced aerodynamic forces for the design of 
active control systens for gust alleviation. 

In theorectical analyses, the change of lift and moment on a wing 
passing through a sharp-edged gust was first calculated for incompressible 
two-dimensional flow by KSrmSn and Sears with simple mathematical formulae 
(Ref. 1). Miles (Ref. 2) extended the calculations to a travelling gust 
field, i.e., sharp-edged gust moving either downstream or upstream rela- 
tive to the airfoil. Drischler and Diederich (Ref. 3) presented results 
for a wide range of wings in both incompressible and compressible flows. 
Meanwhile, the response of an airfoil entering a harmonic gust field was 
first introduced by Sears (Ref. 4). Murrow, et al. (Ref. 5) provided many 
numerical results of lift and moment for finite wings moving through a 
harmonic gust. Giesing, et al. (Ref. 6) also furnished some good sug- 
gestions in computing the oscillatory lift and moment. One notable method, 
called the Doublet Lattice Method (DLM) which was originally developed by 
Albano and Rodden (Ref. 7), was later improved to become a very useful 
tool in unsteady aerodynamics (Refs. P 9). 

For the general harmonic analysis, the atmospheric gust was con- 
sidered as a random set of discrete gusts. Response had been predicted 
most cotnnonly with the assumption that the vertical component of gust 
varied along the flight path, but did not vary along the span. This 
assumption was adopted by most researchers because the associated compu- 
tations were less extensive than these for the more general cases of 



random gusts. It may not be sufficiently accurate for very large aircraft, 
but it should provide useful data for most configurations. 

Besides all these numerical calculations, not much experimental work 
appeared to have been done or to be available. Roberts and Hunt (Refs. 10, 
11) made a series of measurements of transient pressures on a narrow delta 
wing of PR= 1.2 due to vertical gusts, and Patel presented some experi- 
mental results for a couple of delta wings (Ref. 12) and other types of 
wings (Ref. 13) in harmonic gust fields. 

In current aerodynamic research, the vortex flow phenomenon is 
drawing much attention because it offers significant contributions to 
aerodynamic characteristics on low-aspect- ratio wings with sharp or thin 
edges. For these types of configurations, the pressure distributions due 
to the leading-edge vortex separation are drastically different from those 
given by the conventional linear theory. The complex flowfield also makes 
it more difficult to predict aerodynamic forces accurately. Most of the 
aforementioned unsteady aerodynamic methods ar»f based on the linear theory, 
and can not predict the leading-edge vortex effect. Although Atta, et al., 
(Ref. 14) developed an unsteady lifting-surface method with vortex flow 
by using the unsteady vortex- lattice method, no application to various 
gust problems has been presented and the computing cost is expected to be 
very expensive. 

In this report, an unsteady lifting-surface method from reference 15 
will be used to calculate the lift and pitching moment due to sinusoidal 
gusts for several wing planforms. The calculated results are compared with 
other theories for the attached flow and with experimental data for vortex 
flow. In calculating the response of a winy to a sharp-edged gust, Garrick 
(Ref. 16) developed a relation between the oscillatory forces due to flight 
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through continuous harmonic gusts and the indicia! forces due to sharp- 
edged gusts. Instead of Fourier transform taken by most other theories 
in handling this reciprocal relationship, the present method will use Pad6 
approximant to represent the harmonic response and Laplace transform will 
be used to calculate the indicial functions. The problem formulation and 
computed results are presented in the following chapters. 

For many years, military personnel have been continously interested 
in the prediction of the response of an aircraft resulting from a nuclear 
blast wave. Kantian AviDyne (Ref. 17) did a series of experiments to mea- 
sure the blast pressures on a rigid highly sweptback wing at high subsonic 
speeds. McGrew, et al. (Refs. 18, 19) recently used the DLM to develop a 
nuclear blast response computer program for wing-body configurations. Yet 
the method is valid only for the attached flow. The present method, based 
on the unsteady suction analogy, will be used to demonstrate the capabili- 
ty of predicting the nuclear blast response involving the vortex flow. 

Some results in simulating experimental data in reference 17 will be shown. 
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MATHEMATICAL FORMULATION 


2.1 Two-Dimensional Gust Penetration 

Consider that a thin airfoil moving with a velocity U, enters a 
region of atmospheric gust with velocity distribution Wg normal to the 
direction of motion (see Fig. 1). The boundary condition requires that 
the total vertical velocity due to the gust and the vortex sheet simulat- 
ing the airfoil must vanish: 

wg + v/g - 0; for z = 0, -c* < x < c*, (1) 

where w^ is the vertical velocity component of the vortex sheet, x,z are 
the coordinate systems attached to the airfoil and c* Is the semichord 
length of the airfoil. 

To solve a simple harmonic gust problem, the following expression 
is used to specify the sinusoidal gust: 

( 2 ) 

where Wg* Is the amplitude of Wg and w is the radian frequency associated 
with the gust wavelength. Substituting Into Eq. (1) and canceling out 
the time factor e^***^, Eq. (2) leads to 

By way of Introducing the reduced frequency, k = wc*/U, Fq. (3) then be- 
comes 



e-ik(x/c*) 


(4) 


.1 



It Is convenient to divide both sides by U and use the expression x«x/c* 
as the non-dimensional x-coordinate. 

U U 

I.e. , = - Wg e (5) 

where ivg = Wg* / U and = wy*' / U . 

Using Eq. (5), the exact lift distribution for 1ncomp»«ss1ble flow 
can be shown to be (Ref. 20), 

L = 2TTpUc*;^{C(k)[Jo(k) - 1Ji(k)] + Ui(k)} e (6) 

where C(k) Is the Theodorsen's function, Jo(k) and Ji(k) are Bessel func- 
tions of the first kind. For compressible flow, C(k, will be replaced by 
generalized Theodorsen's function** C^Ck) and all other terms In Eq. (6) 
remain the same. 

Furthermore, the lift caused by an arbitrary wg can be calculated 
from Eq. (6). For a sharp-edged gust striking the leading edge of the 
airfoil at t = 0, the boundary condition is 


Wn = 


0 , 


Wo. 


X > Ut - c*. 


X < Ut - c*. 


(7) 


**The generalized Theodorsen's function was generated by Mr. Chung-Hao Hsu 
In his earlier work for calculating lift on wings due to step change In 
angle of attack. 
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Then, Wg can be represented by the Fourier Integral, 




f(w) e dw. 


(8) 


which can be inverted into 


f(w) 




„ -io)t .. 
e dt. 


(9) 


Since Wg = 0 for t < * , Eq. (9) can be shown to be, by following 

=< U 

the Fourier transform of a constant (Ref. 20, p.287). 


f(aj) = 


= f" 


. • _ “lU)t Ux 

Wo e dt 


Woe 


-l 0 )(x+c*) /U 


10) 


(10) 


Substituti«q Eq. (10) back into Eq. (8), the boundary condition for the 
sharp-edged gust Lrcomes 


wo f 

'*3 ' 2 ^ 7 - 


- - T • T> 

CO i(i) 


do) 


» ^ik(s-x-l) 


2^7- Ik 


dk. 


( 11 ) 


where s»Ut/c* is the non-dimensional distance parameter. The airfoil 
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lift due to the harmonic gust (Eq. (5)) is given by Eq. (6). Therefore, 
for the step gust based on Eq. (11), the lift can be calculated as 


L = pUc*v* f ” {C(k)[Jo(k) - iJi(k)] iJi(k)> e 
j, » ik 


ik(s-l) 


dk. (12) 


From Eq. (12), the non-dimensional lift development due to a step 
gust, ^(s), is given by 


L = 2irpUc*Wo\J»(s), 


(13) 


with 




. J_r “ {C(k)lJ«(k)-iJ,(k)l + iJi(k)} 
■ 


ik 


dk. 


(14) 


Let F(k) = C(k)tJo(k)-iJi(k)] + iJi(k), then 


\()(S) = 


■*/ 


ik 


■-/. 




ik 


(15) 


As k is always greater than or equal to zero, Eq. (15) can easily be 
inverted to 


or 


F(k) 


lilli f 

ik 'J 


,|,(s)e'‘*’'^ ds. 


F(k) e"’*' = ikj^ i|;(s) e ds. 


(16) 


» ^ 




^ * 
i ’ 




Letr=ik. Then, 
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= r £ { ip(s) }, 


(17) 


»(s) = L-' { (18) 

r 

where L { i/;(s) } is the Laplace transform of ij^(s). 

Hence, the indicia! lift function can be obtained from the inverse 
Laplace transform involving the amplitude of lift distribution due to a 
sinusoidal gust. 


2.2 Three-Dimensional Gust Penetration 

Now consider a rigid thin wing travellirg at speed U through an 
infinite array of harmonic gusts with vertical velocities Wg, uniformly 
across the wing span. The boundary condition is similar to that for a 
two-dimensional sinusoidal gust: 

w, = . (19) 

where k again is the reduced frequency with a reference chord length c, 
and is a reference point for the gust phase. 

Following reference 21, the oscillatory lift force due to the har- 
monic gust (Eq. (19)) is 

L(t) = pUwgS{L'(k) + iL"(k) }e’’“^. (20) 
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where S is the reference area of the wing, and L'(k), L''(k) are the in- 
phase and out-of-phase components of the dimensionless lift force. 

Next consider a wing entering a sharp-edged gust. Again the 
boundary condition is very similar to Eg. (7), 


' 0, X > Ut + x^g, 

Wg = (21) 

^ Wo, X < Ut + , 

i 6 

where x-jg is the x-coordinate of the leading-edge of the wing. From Eq. 
(20), one can determine the indicial lift function representing dimen- 
sionless lift development due to a step gust, i(<(s), 

L(s) = pUwoS«|/ (s), (22) 

with 

4,(5) . -±f°‘ (.'-'O') * '•■-"(I) (23) 

2tt J _ ik 


Let 


T(k) = L'(k) + if (k). 


(24) 


i|<(s) 



Me^'^^dk. 

ik 


(25) 


With the same procedures from Eqs. (15) - (18), ij<(s) can be obtained from 
the inverse Laplace transform involving the dimensionless lift distribu- 
tion due to a harmonic gust, 

,|,(s) = {JM }. (26) 

r 
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2.3 Pad^ Approximant 

To be able to calculate the inverse Laplace transform indicated in 
Eqs. (18) and (26), it is convenient to express T(k) and F(k)e as closed- 
form functions of r. In any lifting-surface computation, T(k) or F(k) are 
calculated only at a finite number of k's. These values can then be in- 
terpolated by Pad^ approximant as suggested by Vepa (Ref. 22). Following 
Vepa, an [ N,N ] sequence of Pad^ approximant to approximate T(k) in three- 
dimensional case ano F(k)e"’^ in two-dimensional gust can be written as, 

T(k) = 1 - ( N,N ], F(k)e"'" = 1 - [ N,N 1, (27) 

( N,N I = A.rW t A,r«-1 * t Aw-lr _ (28) 

r^ + A|^r^~^ + + 

where A-j are the coefficients of Pad^ approximant. 

For N = 3 and use T(k) as an example. 


T(k) = 1 - A,rH A,r^-> A;r ^ , 2 ,,) 

r^+ Aar^+ A4r + A5 

and for N = 2, 

T(k) = 1 - . (29b) 

r"+ A 2 r + As 

Let 

1 - T(k) = V(k) + iW(k). (30) 

From Eq. (30), it follows that for N = 3, 
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V(k) + iW(k)= - ^ \ (31) 

(ik)^+ A3(ik)*+ A4{ik) + A5 , 

Eq. (31) can be expanded and separated Into real and imaginary 
!>arts respectively to give, 
real part; 

Aik^- AjVk- A 4 Wk + AsV + Wk^= 0, (32a) 

imaginary part: 

Aok^- Aak- AsWk + A^Wk + AsW- Vk^= 0. (32b) 

Normally, more values of k's are chosen to calculate T(k) than the 
number of unknowns A^ in Eq. (32a) and (32b). Therefore, a least square 
technique to be described next is used to determine Aj's. 


2.4 Least Square Technique 

The Ip st square principle is based on the requirements that A^'s 
are determined by minimizing the sum of squares of errors: 

M M 

Sum= ^ (L.H.S. of Eq. (32a))^+ ^ (L.H.S. of Eq. (32b) )^. (33) 
1 1 

where M is "e number of k's of which T(k) is calculated by any existing 
1 i f ti nC'Surface theory. 

At a minimum, all the partial derivatives with respect to A/s, 
;ch a"> 3Sum/3Ao, 3Sum/3Ai, , 3Sum/3A^, must vanish. These condi- 

tions result in i+1 equations for i+1 unknown coefficients A^-'s. Thus, 
A^'s c'.n je solved from Eq. (33) with M different lift values at corre- 
sponding k's. 
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2.5 Indicia! Lift Function 

After all coefficients of the Pad6 approximant are determined, Eq. 
(29a) can be rewritten as. 


T^k) = 1 - r AprH Air -»• Az 

r’+ Asr^+ A4r + As 

By partial fraction method, Eq. (34) leads to 

T(k) _ J_ f Bj 

r r - Bi ’ 


(34) 


(35) 


where is the ith root of the polynomial in the denomenator of Eq. (34) 
and B-j is the corresponding coefficient of partial fraction in Eq. (35). 

Based on Eq. (26), the indicial lift function i/;(s) can be obtained 
by applying the inverse Laplace transform to Eq. (35), 


♦(s) = 1 - E 
i=l 


(36) 


2.6 Nuclear Blast Response 

Calculation of the lift development of a thin wing encountering a 
nuclear blast wave will follow the same way as in computing the sharp- 
edged gust response. No major change has to be made except the boundary 
condition. The shock wave induced by the nuclear blast is assumed travel- 
ing at sonic speed. Thus, an impulse function 6(s) is used instead of the 
unit step function in the gust response condition (Ref. 24), 

X - Xi- 

Wg = Wo6(t = ). (37) 
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where U is the magnitude of the vector sum of the shock wave velocity and 
the freestream velocity. Eq. (9) now becomes, 


f(w) 


■/.: 


Wo6(t - 


X - X 


li.) e -’“‘dt 


-iw(- 


X - X 


le 


) 


(38) 


= Woe ' U 

Comparing Eq. (38) with Eq. (10), the main difference is the factor "ioj" 
in the denomenator of Eq. (10). This also follows from the fact that the 
impulse is equal to the time derivative of the step input. Hence, 

P{h(t)} = i(uF{A(t)} , (39) 

where h(t) is the unit impulse, A(t) is the step function and F{ } is 
the Fourier transform. 

Rewriting Eq. (25) for the impulse response, it follows that. 


h(s) 




T(k)e ^'^^dk. 


(40) 


or. 


h(s) = L~' { T(k) } 


(41) 


The Pad^ approximant used for the indicial response analysis re- 
mains applicable, except the inverse transform calculation has changed. 
For example, N = 2 Pad^ approximant can be written as. 


T(k) .. 1 - 

r^+ A^r + As 
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2 B. 

= 1 - Ao - £ ~Q - . (42) 

i=l 


From Eq. (41), the impulse lift function h(s) can obtained. 


h(s) = {T(k) } 


= (1 - A )6(s) - 


I Bi 
i=l 


(43) 
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3. NUMERICAL RESULTS AND DISCUSSIONS 


3.1. Gust Penetration 

3.1.1. Sinusoidal gust problem 

The lift development due to a harmonic gust is calculated by the 
computer program based upon the unsteady quasi-vortex-lattice method de- 
veloped by Lan (Ref. 23). 

For a thin airfoil, 30 vortex elements are used in the computation. 
The steady state two-dimensional lift values are simply calculated by the 

equation Ct = 2ti/6, where 3 =/l-M^, M is the frees tream Mach number. The 

*0 

computed results are compared with Sears' for incompressible flow and with 
Graham's for several different Mach numbers at reduced frequency 2. Both 
comparisons, showing good agreement, are tabulated in Tables I and II. 

The three-dimensional unsteady aerodynamic program of reference 15 
is then revised to account for the gust response. The present attached 
flow results of lift and moment for a delta wing of 75° sweep at M = 0.4 
are compared in Figs. 2 and 3 with those calculated by a kernel function 
method (Ref. 5). 

Results of calculation will also be compared with experimental data. 
In references 12 and 13, a gust tunnel which could generate a sinusoidal 
vertical gust was used to measure the oscillatory lift and moment on two 
delta wings of AR=1 and 2, and several other commonly used wing planforms. 
The tests were performed for all wings at two mean freestream velocities 
of 12.43 and 20.00 m/s. However, the oscillatory gust waves convected 
downstream with a velocity of 0.61 of the mean freestream velocity. This 
would indeed influence the gust wavelength, i.e. , the frequency parameter. 
Thus, in the present calculations based on gust moving with the freestream 
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velocity, an effective frequency (k' » k/0.61) , as suggested by Patel, will 
be used in the following comparisons. 

In Fig. 4, test data for a rectangular wing of 4R = 6 (Ref. 13) are 
compared with two sets of theoretical results. It is seen that the pre- 
sent theory agrees well with Graham's (Ref. 24) in the predicted in-phase 
component of the oscillatory lift. The phase lag is underpredicted and 
both theories overpredict the lift amplitude. 

To demonstrate the vortex flow effect, a delta wing of /R=l is used 
to compare with Patel's data (Ref. 12) at a = 0 and 12 degrees in Fig. 5. 
The present theory predicts well the lift amplitude at both angles of at- 
tack. At a=12°, the phase lag is also adequately predicted by the pre- 
sent theory. On the other hand, the phase lag at a=0° is not accurately 
predicted. In seeking the reasons for these deviations, several points 
should be noted: 

(1) At a=0 degree, the attached flow prevails. It is of interest 
to compare the present results with the doublet-lattice method (Refs. 8,9) 
for this delta wing of ^=1 . In Fig. 6, agreement between the present 
results and doublet-lattice method's is excellent in lift amplitude and 
phase angles. 

(2) As depicted in reference 12, force measurements were made re- 
lative to the undisturbed frees tream gust at the root 2/3rd chord point. 
Another point at the gust tunnel exit was also used as a reference. It 
is not known whether the conditions at the exit point were disturbed once 
the test model was placed in the tunnel. 

(3) Also, Patel indicated in reference 12 from test data that the 
incremental lift due to vortex lift contribution was important in magni- 
tude only with no measurable contribution to phase angles. However, the 



present theory shows that this is approximately true only with respect to 
some reference points li.e., Xp at Eq. (19)]. This is illustrated in Fig. 
7 with the root midchord point and wing apex as reference points. There- 
fore, the results by the present theory very much depend on the precise 
location of the reference about which the phase angle is calculated. 

At any rate, reasons for the discrepancy in the predicted phase 
angles at o=0 degree for the delta wing of AR = 1 are still unknown at 
the present time. 

3.1.2. Fade approximant 

Two different sequences of Padd approximant are constructed here 
to fit various freestream conditions. For airfoils in compressible flows, 
Pad6 approximant with N = 3 in Eq. (27) (called Pad^ A6) is used. On the 
other hand, Padd approximant with N=2 (called Fade A3) is used for air- 
foils in incompressible flow and also for three-dimensional conditions. 
These choices are made through numerical correlations. The corresponding 
matrices for deciding the unknown coefficients A-j in Eq. (28) are present- 
ed in Appendix A. 

The good agreement between the calculated results from the unsteady 
QVLM and by Fade approximant for a thin airfoil in the harmonic gust is 
shown in Fig. 8 for incompressible flow. Fig. 9 for Mach 0.5 and Fig. 10 
for Mach 0.7. The oscillatory lift for a delta wing of AR = 1.2 to be used 
for calculating step gust response is shown in Figs. 11 and 12 for upward 
and downward gust respectively. 

Most other theories use Fourier transform to handle the reciprocal 
relationship between the oscillatory and indicial lift forces. So they 
need harmonic lift values at a large number of different frequencies. 
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But the number of k required by the present method is only 7 (including 
k = 0), and the indicial lift results, which will be seen in the next sec- 
tion, still show good accuracy. 

3.1.3. I.dicial lift function 

The indicial lift results performed by using the inverse Laplace 
transform described in the preceding chapter will be presented. First, 
for a thin airfoil passing through step gusts in incompressible and com- 
pressible flows. Fig. 13 shows plots for Kiissner function, and Tigs, 14, 

15 exhibit the indicial lift at Mach 0.5 and 0.7. Oood results have 
been expected because of the accurate approximant shown in last section. 
The exact solutions are calculated through the data in Table 6-2 of refer- 
ence 20. From Fig. 16, it is apparent that the compressibility effect 
decreases the rate of lift build-up in two-dimensional flow. 

Second, in thin wing gust penetration, for lack of experimental 
force data to make a direct comparison, some pressure data (Ref. 11) will 
be employed to compare the trend produced by calculated total lift. 

The configuration of the delta wing used in reference 11 is shown 
in Fig. 17(a). The model wing was carried along a straight railway track 
on a rocket-propelled sledge, through the efflux from an open- jet wind 
tunnel blowing across the track. The velocity of the sledge was 180 ft/ 
sec., the tunnel efflux velocity was 47 ft/sec., and the model was at zero 
angle of attack. For measurements of transient pressure, four transducers 
were positioned at locations being 0.3, 0.4, 0.5, 0.6 root chords aft of 
the apex, and along a line at 75% semispan [see Fig. 17(a)]. In accord- 
ance with figures in reference 11, the indicial lift is plotted against 
local distance parameters Xp/b^ in comparing with pressure values through 
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transducers A, B and C, where and are defined in Fig. 17(b). 

Fig. 18 shows the similar trend between measured transient pressure 
and the calculated indicial lift for both upward and downward sharp-edged 
gust. It is seen that the development of a gust-induced gain of lift is 
very gradual. On the contrary, the gust- induced loss of lift occurs 
relatively instantly. In the calculations, the vortex lag descrii>**d in 
reference 15 is assumed to be present if the lift is increasing and th?re 
is no vortex lag if the lift is decreasing. This assumption appe; o be 
reasonably accurate for the leading-edge vortex flow. This can ai ' 
seen from Fig. 19 which illustrates the comparisons among the vortex flow 
and potential results. There is some significant difference for the up- 
ward gust while the trend is quite close in the downward gust condition. 

3.2. Nuclear Blast Response 

The unsteady aerodynamic program of reference 15 is again used to 
calculate the nuclear blast responses of aircraft flying at high subsonic 
speeds. 

Reference 17 is the only obtainable test data which can be used to 
check the leading-edge vortex separation effect, predicted by the present 
method, on aircraft nuclear blast response. The sideview of a thin wing 
intercepted by nuclear blast waves is shown in Fig. 20. Fig. 21 shows the 
test model which consists of a swept wing of 67° leading-edge sweptback 
angle with a nose and partial fuselage section. In the test, the model 
was mounted on a high speed dual rail rocket sled at an initial angle of 
attack of 3.2° . The sled, travelling at Mach 0.76, was intercepted pro- 
gressively by blast waves from sequential detonation of charges of TNT 
with the blast intercept angle 9 = 20° (see Fig. 20). Twenty pairs of 


19 


pressure transducers were Installed on the left wing half to measure the 
blast- Induced pressures. 

For the purpose of correlation with the test model, a semiwing used 
by the present method Is Illustrated In Fig. 22. For such a configuration, 
the concept of augmented lift Is Included In the present calculation. The 
definition for the characteristic length Is adopted from reference 26. With 
negative augmented vortex lift, the vortex lift effect may not be as strong 
as expected even the leading-edge sweptback angle is 67° for the test model. 

Fig. 2'', reproduced from reference 17, shows the pressure variation 
measured by transducer 13 which was positioned at half semi span and along 
quarter chord line. Like DLM, the present theory underpredict the blast- 
induced incremental pressure loadings because the nonlinear vortex effect 
is not Included in the calculation of pressure differential. The present 
theory, being based on unsteady suction analogy, can only demonstrate the 
vortex effect by the variation of total lift or moment. Fig. 24 shows the 
comparison of vortex flow with attached flow for impulse lift. It Is evi- 
dent that the vortex lag decreases the rate of lift decay after the blast 
Intercept. 

There are several factors which must be mentioned in connection with 
Figs. 23 and 24: 

(1) The blast amplitude is determined by peak material velocity 
(gust) behind the shock front and Is assumed invariant throughout the 
present calculation. 

(2) In the experiment, because of different locations of transducers, 
the blast Intercept time (time at the shock arrival) Is different for each 
transducer. In the present calculation, the coordinate origin Is set at 
the apex of the wing, so that the Intercept Is exactly at S“0 as 
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shown in Fig. 24. However, to correlate with the test data, a shift has 
been made In Fig. 23. 

(3) The steady state lift value is used for all time parameters 
less than zero since there is no incremental lift or pressure differential 
values for s < 0 in the present analytical results. 

For a more significant vortex lift effect, the impulse lift on a 
delta wing of 4R = 2 at Mach 0.5 intercepted by a blast wave with the same 
intercept angle as that in Fig. 20 is considered. The initial angle of 
attack is assumed to be 15** and the peak blast-induced angle of attack 
is assumed to be an additional 15*. The attached-flow oscillatory lift 
can be well represented by a Pad6 approximant as indicated in Fig. 25. 
However, with the vortex lift effect included, the Pad6 approximant 
fails to approximate the calculated results accurately as shown in Fig. 

26. Because of this reason, the calculation of impulse lift for this 
delta wing was no*" successful; and hence, the results are not presented 
f .re. 
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4. CONCLUSIONS AND RECOMMENDATIONS 

Art unsteady lifting-surface computer program based on quasi-vortex- 
lattice method along with leading-edge suction analogy has been developed 
to estimate the oscillatory air forces on wings of general planforms in 
gust flow at any frequency. Pad€ approximant and Laplace transform have 
made it practical to convert the oscillatory air forces to indicial air 
forces . 

Both the experimental data and other theoretical results are used 
to check the accuracy of the present calculations in attached flow and 
with vortex lift effect. It is shown that the present method can accu- 
rately predict the oscillatory and indicial lift on wings in different 
gust fields. Also, the phenomenon of the gust-induced gain of lift being 
very gradual and the gust-induced loss of lift occurring relatively ab- 
ruptly can be explained by the presence or absence of vortex lag effect. 

The present program is extended to acco'jnt for the nuclear blast 
response as well. Though there is no lift data available at this time, 
the trend for the vortex lag is clearly seen from the comparison between 
the attached flow and vortex results predicted by the present method. 

The following points should be noted to improve the efficiency 
and capability of the present method. 

(1) In the present calculations, 72 vortex elements were used for 
half wingspan - 6 in the chordwise direction and 12 in the spanwise di- 
rection. It is recommended that 40 elements (e.g., 4 in the chordwise 
direction, 10 in the spanwise direction) could be used in lieu of 72 for 
small reduced frequencies; thus, the size of the aerodynamic influence 
coefficient matrice could be greatly reduced. 
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(2) The present method can only deal with nuclear blast locations 
being on the plane of symmetry, i.e. , the X-Y plane. For more general 
blast orientations, the present method should be extended to calculate the 
aircraft wing response of the blast waves coming from any arbitrary direc- 
tion, since an asymnetric condition can always be treated as a combination 
of a symmetric and an anti symmetrical conditions. 

(3) Among the existing approximation schemes for the unsteady aero- 
dynamic loads, Karpel's approximant (Ref. 28) provides better accuracy than 
Fade's. It is suggested the present method should include more than one 
approximant to meet different gust and blast conditions. 
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Table I. Oscillatory Lift* on an Airfoil in Sinusoidal 
Gusts at M = 0. 



Presen 

t Method 

Sears 

j ** 

k 


^li/^io 

^ir/^io 

^li^^lo 

0.02 

0.96289 

-0.07484 

0.963 

-0.075 

0.04 

0.92403 

-0.11449 

0.924 

-0.114 

0.06 

0.88697 

-0.13923 

0.887 

-0.139 

0.08 

0.85266 

-0.15457 

0.852 

-0.154 

0.10 

0.82126 

-0.16348 

0.821 

-0.164 

0.20 

0.70156 

-0.15964 

0.702 

-0.160 

0.40 

0.56789 

-0.08494 

0.568 

-0.085 

0.60 

0.48837 

-0.00490 

0.488 

-0.005 

1.00 

0.36865 

0.12594 

0.369 

0.126 

2.00 

0.08158 

0.26796 

0.082 

0.268 

1 ^ and 

C.|^ are the 

in-phase and out-of-phase 

components 


**Sea*"s' results are copied from reference 27. 
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Table II. Oscillatory Lift on an Airfoil in Sinusoidal 
Gust at k = 2.0. 


Present Method Graham's 


M 

^ir 

^li 

^Ir 

^li 

0.0 

0.5126 

1.6836 

0.5125 

1.6837 

0.2 

0.6925 

1.7473 

0.6955 

1.7441 

0.4 

1.1535 

1.5960 

1.1511 

1.5953 

0.5 

1.2773 

1.3139 

1.2774 

1.3138 

0.6 

1.2U8 

1.1104 

1.2139 

1.1102 

0.7 

1.2192 

1.0113 

1.2152 

1.0085 

0.8 

1.1881 

0.8921 

1.1838 

0.8897 

0.9 

1.1516 

0.7986 

1.1486 

0.7961 


The results correspond to v^=l in Eq. (4). 

•kit 

Graham's results are taken from reference 23. 
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APPENDIX A Coefficient Matrices of Pad6 Approximation with Least Square 
Techniques 
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Fig. 3 Generalized oscillatory pitching moment for a 75® delta 
wing due to harmonic gusts with Wg = 1 at M»0,4; pitching axis 
at root midchord point. 
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Fig. 6 Incremental lift for a delta wing of /R = 1 in sinusoidal 
gusts for a *0® with w = 0.0314 and reference point at root 2/3rd 
chord point. 
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Fig. 7 Variations of phase angle due to changes of angle of attack 
and reference points for a delta wing of = 1 in sinusoidal gusts. 
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iscillatory lift for a thin airfoil due to hannoni 






Fig. 10 Oscillatory lift for a thin airfoil due to hannonic 
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Oscillatory lift for a delta wing of /JR = 
al gusts for a=0° with Wq = 0.261 at M=0. 
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dilatory lift for a delta wing of y<R = 1.2 in 
gusts for a= 14.6° with ^ = -0.261 at M= 0.16. 
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Fig. 13 Indicia! lift function from a sharp-edged gust of a thin airfoil at M= 





Indicia! lift functi 
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Fig. 17(a) Gereral configuration and transducer positions of the 
test iBodel In reference 11. 
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Fig. 18 Trend of indicia! lift function from vertical gusts of a delta wing of 




Fig. 18 Continued. 
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Fig. 18 Continued. 
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ft function for a delta wing of /IR = 1.2 entering into vertical gusts at M=0.16. 
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Fig, 20 Thin wing intercepted by a blast wave moving with sonic speed. 
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Fig. 21 General configuration of the test model used for nuclear blast | 

response (taken from reference 17). j 
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Fig. 22 Simplified wing planform us 
for nuclear blast response analysis. 
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Fig. 23 Pressure variation of the test wing due to nuclear blast waves with 9*20® for a* 3.2' 
at N»0.76 (reproduced from reference 17). 
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Fig. 24 Impulse lift function of the test wing due tc nuclear blast waves with 0 = 20° for 
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